Quantum superconductor-metal transition in a proximity array 
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A theory of the zero-temperature superconductor-metal transition is developed for an array of 
superconductive islands (of size d) coupled via a disordered two-dimensional conductor with the 
dimensionless conductance g = %/e 2 Ra S> 1. At T = macroscopically superconductive state of 
the array with lattice spacing h 3> d is destroyed at g < g c ~ 0.1 ln 2 (6/d). At high temperatures the 
normal-state resistance between neighboring islands at b — b c is much smaller than Rq — h/4e 2 . 

PACS numbers: 74.40. +k 71.30.+h 74.50. +r 



In two-dimensional (2D) systems, qualitative argu- 
ments based on duality between Cooper pairs and vor- 
tices lead to the prediction |3J] that the superconductor- 
insulator transition happens at the universal quantum 
value Rq = h/4e 2 of the resistance per square R n . Al- 
though a number of experiments (cf. [0 and references 
therein) seem to be in agreement with this prediction, 
other data demonstrate strong deviations from it HI). 
Phenomenological picture of duality is not able to predict 
the system's parameters (e. g., the value of the normal- 
state resistance) leading to the quantum critical point: 
a microscopic theory is needed to find it. Competi- 
tion between Josephson coupling Ej and charging en- 
ergy Eq is known |t],^| to be the driving mechanism of 
zero-temperature phase transitions between the super- 
conductive and insulating states in artificial arrays |^],||, 
films jlC],[llJ and bulk materials |l2|. In such systems 
there are no free electrons at very low temperatures due 
to Cooper pairing, but pairs may become localized due 
to Coulomb repulsion. This is the "bosonic" mechanism 
of superconductivity suppression. 

Homogeneously disordered superconductive films 
JOHHI P resen t another group of systems where quan- 
tum fluctuations lead to destruction of superconduc- 
tivity. The theory of T c suppression in such films was 
developed in Refs. JlJ|. The qualitative idea behind 
this theory is that disorder-enhanced Coulomb repul- 
sion leads to the decrease of Cooper attraction and 
thus to the decrease of T c . The superconductive tran- 
sition temperature vanishes [p^JT^ I when the dimen- 
sionless film conductance g — h/e 2 Ra decreases down 
to g Fhl = (27r)- 2 ln 2 (l/T c0 T tl .), where T c0 is the BCS 
transition temperature and r tr is the elastic scattering 
time. This second ( "fermionic" ) mechanism of supercon- 
ductivity suppression is clearly different from the first 
one since its basic feature lies in disappearance of 
Cooper pairs altogether. Experimental data supporting 
the fermionic mechanism are reviewed in Ref. paL A 
drawback of this theory is that it neglects quantum fluc- 
tuations of the bosonic field (i. e., it can be considered 
as a kind of the BCS theory with the renormalized at- 



traction constant). For a phenomenological comparison 
of the bosonic and fermionic mechanisms, cf. Ref. j^ j. 

In this paper we study a model for quantum break- 
down of superconductivity, which lies in between the two 
limiting cases discussed above. We consider an array of 
small superconductive (SC) islands (of radius d each) in 
contact with a thin film of dirty normal (N) conductor 
with the dimensionless conductance j> 1. The distance 
between neighboring islands is 6 > d (more precisely, 
b~ 2 will be the concentration of islands). Resistance Rt 
of the interface between each island and the film is low: 
Gt = h/e 2 RT ^S> 1. Islands are thick enough, to prevent 
suppression of superconductivity inside them. The cor- 
responding condition for the superconductive gap reads 
Asc 3> Gt/vVi, where Vi is the island's volume and v 
is the density of states. The above inequality ensures 
that the lifetime of electrons in the SC island (which is 
finite due to the tunneling across the interface) is much 
longer than %/ A$c- We assume also that G\ ^> Ang; the 
meaning of this condition will be explained below. We 
will show that macroscopic superconductivity in such a 
system at T = becomes unstable with respect to quan- 
tum fluctuations at g less than 



9c = Gc 



7T d 



(1) 



where Q c *~ 1 will be determined below, and d ~ d. 
Eq. (p]) presents our main result (obtained within log- 
arithmic accuracy), which shows that the critical sheet 
resistance Ra c = ^/e 2 3c is much less than the quantum 
resistance Rq, provided ln(b/d) > 3. Moreover, the same 
is valid for the normal-state resistance between neigh- 
boring islands R n — (Rn c /^) ln(b/d) ~ 2>Rq / \n{b / d) . 
This result is at odds with usual arguments based on the 
model of resistively shunted Josephson junctions : in 
that model superconductive behavior of a single junc- 
tion is preserved at T — as long as it has resistance 
R < Rq. The physical reason for this discrepancy is 
that we account for the discrete nature of charge trans- 
port between the SC islands, which was neglected in the 
Caldeira-Leggett model adopted in p7|. 
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We will follow an idea presented in Ref. p8[ , where a 
simplified version of the considered model was analyzed 
(cf. also f0|). Namely, we make use of the long-range 
nature of the Josephson coupling Jy between SC islands 
due to the proximity effect in the film, which scales as 
Jij oc r~j 2 exp[— c(rij / Lt)], where Lx = \JKD/T is the 
thermal coherence length (cf. Refs. jl8|,^0|), c ~ 1, and D 
is the diffusion constant in the film. At low temperatures 
the interaction radius Lx diverges indicating that the po- 
sition of the quantum phase transition can be found in 
the mean-field approximation (MFA) analogous to the 
one developed in Ref. |23| . Within the MFA, macroscopic 
superconductive coherence sets in at 



Lj(T)C{T) > 1, 



where 



J{T) = V J tj , C(T) = [ 1/T drCo(r), 

Jo 



(2) 



(3) 



and Co(t) — (cos((9(0) — 9(t))} is the single-island auto- 
correlation function of the order parameter phase. Thus, 
to find the superconductor-metal transition line we need 
to determine the functions J(T) and C(T) at low tem- 
peratures. We will discuss the T — limit mainly, and 
present the results for low T > at the very end of the 
paper. 

Calculation of J{Q) is rather straightforward. The re- 
sult depends on the relation between the dimensionless 
Cooper repulsion constant A„ defined at the energy scale 
LUd = D/d 2 and the ratio Gx/^g which is a measure of 
the interface transparency. At low interface transparen- 
cies, Gr/47rg <C A„, it is possible to consider the lowest 
order process over Gt, and one obtains [[L8[ 



J(0) 



1 



b 2 ln(b/d) ' 



(4) 



This result is similar to the one obtained in j2(J for a one- 
dimensional S-N-S geometry: in both cases J{0) formally 
diverges as the Cooper repulsion constant goes to zero. 
However, the decrease of A„ at a fixed Gt drives the sys- 
tem out of the low-transparency limit. At A„ < Gx/^g 
full nonlinear treatment of the proximity problem is nec- 
essary, accounting multiple Andreev reflections from SC 
islands. At the same time, the effect of A n upon JT(O) can 
be neglected in this limit (for 2D geometry studied here) . 
To find the coupling energy J(r) between two islands sit- 
uated at the distance r ^> d, one has to solve the Usadel 



equation |21| in the film. Actually, it is more convenient 
to calculate the maximum superconducting current be- 
tween the islands, I c (r) = (2e/%)J(r). The value of the 
supercurrent can be determined as the integral over the 
current density, taken along the straight line which lies 
in the middle between two islands, and is perpendicular 
to the vector r connecting them. If the distance between 
the islands is very large, m(r/eQ ^ 1, anomalous part of 
the Green function is small along the above-mentioned 
line, that allows 1 to linearize the Usadel equation at rel- 
evant energies E ~ TiD/r 2 and find Icif)- Thereby one 
obtains J(r) = (n 3 /A)gD/r 2 \n 2 (r/d) and, finally, 



J(0) 



gD 



2 b 2 \n(b/d)' 



(5) 



The key point in the discussion of the T — transition 
is to determine C(T — > 0). We will see that C(0) depends 
exponentially on the film conductance g. If islands do 
not have ohmic contacts with the film (coupling via ca- 
pacitance Cj only) then C(0) = h/Ec — hCj/2e 2 . In our 
case h/C(0) plays the role of an effective charging energy 
Eq of an island that survives in spite of good conduc- 
tance around. To make ideas transparent, we first discuss 
a simplified model [jl8| with sufficiently strong Cooper- 
channel repulsion in the film, A„ 3> Gx/^g- Then dy- 
namics of the phase 9(t) of a single SC island can be 
described by a simple imaginary-time action, 



(I 



(t _ t0 2 ■ (6) 



Here Ga — Gx/4irg\ n is the Andreev subgap conduc- 
tance (normalized to e 2 /%) in the limit of weak proximity 
effect, valid under the condition A„ 3> Gx/^g |l8|j24[| 2 . 
Expression m) is valid at low energies, u> < Wd e_1 j 



while at higher energies, o^e 



-l/A* 



< lu < ujd, one has 



Ga(w) = (G T /4:TTg)hi(ui c i/^)- Thus, the Andreev con- 
ductance is large compared to Ae 2 /h for all frequencies 
below iOd if the condition G\ ^> Aug is fulfilled. At lower 
values of Gx, N-S transport across the interface is sup- 
pressed by the usual Coulomb blockade effect governed 
by the junction's charging energy 2e 2 /Cj. We are not 
interested in this trivial effect, and will assume the con- 
dition G\ ^> 4wg to be valid. 

For large Ga one can start from the Gaussian approx- 
imation for So[0{t)]. Then the Fourier-transformed cor- 
relator of phase fluctuations (|# w | 2 )o = 4/\u>\Ga, and, 
hence, C a (t) = e -K(0«-0(o)) 2 >o oc t^^ GA . At G A > 
4/-7T, C(T — > 0) diverges that seems to indicate that at 
large Ga superconductivity is always stable at T — 0, in 



1 Similar calculation can be found in Appendix C of Ref. |22] 

2 The result for Ga given in JI3] is overestimated by the fac- 
tor of 2 due to a numerical mistake. 
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agreement with p7jj. The crucial point is to note that 
the employed Gaussian approximation breaks down at a 
finite time scale t* , due to downscale renormalization of 
G a- This renormalization is caused by the periodicity of 
the action So[9] as a functional of 8(t), that is, in physical 
terms, by the charge quantization. This problem is anal- 
ogous to the one studied by Kosterlitz | p5[ . Translating 
his results to the present case, one gets the renormal- 
ization group (RG) equation dGA(()/d( = — 4/7r, with 
Q = lna^i- This equation is to be solved with the initial 
condition Ga(0) = Ga- As a result, at the time scale 
t* ~ w7 1 e 7rG - 4 / 4 the renormalized Andreev conductance 
Ga [t* ) decays down to the value of order unity jl8| . At 
longer time scales Co(i) decays approximately as t~ 2 , so 
the integral C(0) ~ t* ~ Lu^ 1 e 7rGA ^. Taking into ac- 
count that J ~ b~ 2 , and using Eq. (||), one obtains |ll| 
the critical distance between islands b c ~ de 7rG - 4 / 8 . 

However, this result is valid under the condition X n 3> 
Gx/^g which is difficult to realize simultaneously with 
the inequality G\ ^S> Ang needed to avoid trivial Coulomb 
blockade effects. Indeed, at energies E <C h/T tI , Cooper 
interaction constant X(E) is determined by the RG equa- 
tion Q which we present in a simplified form p4j | valid 
for hi(h/ Er tr ) <C g when renormalization of g can be 
neglected: 



dX - \ 2 4- \ 2 



(7) 



and X(( = 0) = A„. The fixed point solution of Eq. (]?]), 
A = X g , is too small to fulfil both the above inequali- 
ties together. Therefore, typically the approximation of 
single-parameter RG for Ga is not valid, and we should 
reconsider the problem of the subgap N-S conductance 
in the presence of three different effects acting simulta- 
neously: 

(i) disorder-enhanced multiple Andreev reflections |2(| 
which increase Ga(^) when the coherence length 
\J~DJui grows; 

(ii) Cooper-channel repulsion A which reduces 
G A PWPl; 



(iii) quantum fluctuations of the phase 9(t) which de- 
stroy coherence between Andreev reflections and 
suppress Ga(() at long time scales. 

To treat all these effects together, we employ the func- 
tional RG method for the proximity-effect action in the 
Keldysh form @. 

Like in the simplified model jl8| discussed above, the 
constant C(0) is determined (with exponential accuracy) 
by the value of time t* when Ga(( — hic^i*) becomes 
of the order of 1, since at longer times Co(t) decays fast. 
However, the equation for Ga(C) is much more compli- 
cated now as it includes an infinite set of parameters. To 
derive the corresponding RG equations, we start from the 



Keldysh action for a SC island in contact with a disor- 
dered metal, derived in Ref. pi}] . It can be represented 
as a sum S = S'buik + Abound of the bulk and boundary 
[the last term in Eq. (0)] contributions: 
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D(VQ) 2 + Ai(ir z d t + <f> + A)Q 



2v . 



■Tr ^V' 1 <f> + ±f Tr A+a x A 
A 



ittGt 



TrQsQ. (8) 



The bulk action, S'buik, is a functional of three fluctuat- 
ing fields: the matter field Q(r,t,t') in the film (its av- 
erage value gives the time-domain representation of the 
electron Green function G(r,r') at r = r'), the electro- 
magnetic potential 0(r, t), and the order parameter field 
A(r, t) used to decouple the quartic interaction vertex in 
the Cooper channel. Q(r, t, t') is a matrix in the time do- 
main, and in the direct 4 x 4-dimensional product K®N 
of the Keldysh and Nambu-Gor'kov spaces. Pauli ma- 
trices in the K and N spaces are denoted by cr, and Tj, 
respectively. The field Q satisfies a nonlinear constraint 
Q 2 = 1 and can be parametrized as Q = e~ w / 2 Ae w / 2 
with {IT, A} = 0, where A = Aot z is the metallic saddle 
point and 



Ao(e) 



1 2F(e) 
-1 



(9) 



A' 



The matrix F(e) — To/(e) + r z /i(e) has the meaning 
of a generalized distribution function. The object cf> = 
(4>ii 02 ) T is a vector in the Keldysh space, with <fii, </> 2 be- 
ing the classical and quantum components of the ^-field. 
4> is a shorthand notation for the matrix <f) — 4>i(Jo+(j)2<J x . 
Similarly, A = (Ai, A2) T , and A stands for a 4 x 4 ma- 
trix A = [t + Ai - r_A^](To + [t+A 2 - t^A* 2 ]cj x , where 
t± = (t x ± it v )/2. In terms of the tr-model action (||), 
diffuson and Cooperon collective modes of the electron 
system are described as slow fluctuations of the Q-matrix 
over the manifold Q — 1. The last (boundary) term in 
Eq. (||) describes an elementary tunneling process be- 
tween the SC island and the N metal. The matrix Qs 
describes the state of the SC island. At the low-energy 
scales e <C \Asc\ it is expressed via the phase 6{t): 



Qs = -ir+e l 



(10) 



where 9 = 8i<7q + 92<J X . 

The action (||) contains a fluctuating scalar potential 
field (j> accounting the direct Coulomb interaction in the 
density-density channel. Major effects of this interac- 
tion are: (i) local electroneutrality of electron liquid at 
low frequencies, and (ii) zero-bias anomaly in the tun- 
neling DOS pq] . Both effects can be taken care of by 
means of a special gauge transformation E9]: Q t t> — > 



a iK{t)r z 



-iK(t')Tz 



and <f>(t) -> + d t K(t). The 
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"Coulomb phase" K(t) is a linear functional of <f>(t) cho- 
sen in a way to cancel all terms which are linear both in 
4>{t) and in the matrix field W(r, t,t'). Then the effect 
(i) is contained in the tree level of the transformed ef- 
fective action, whereas (ii) comes from the simplest loop 
correction j2^j24| . After the above gauge transformation 
the phase 9j(t) of the j-th island enters the action in the 
combination 9 j{t) — 2K(t,rj) only. Now the important 

point of our discussion comes about: the phases 9j(t) 
of each island are not fixed by any external source, and 
should be integrated out. Thus the shift of integration 
variable 9 At) — > 9j(t) — 2K(t,rj) eliminates K(t,rj) 
from the action, together with both effects (i) and (ii). 
In other terms, the present problem of unconstrained 
phase 9 (t) fluctuations can be treated as if it would be 
no Coulomb interaction, since Gaussian terms in the ac- 
tion containing electric field are decoupled from the rede- 
fined 9 (t) variable. It is thus legitimate to neglect elec- 
troneutrality and calculate frequency-dependent subgap 
conductance Gyi(u) as if the outer normal contact would 
be placed at the distance = \/D/lo from the SC is- 
land. We emphasize that the same would be wrong for a 
usual problem of N-S conductance between contacts with 
fixed voltages, where the full size of the N film, L ^> 
does enter the result, adding the term (Rn/lir) ln(L/i? w ) 
into the resistance, cf. Sec. VI of Ref. p2]| . 

Next we use the RG method to integrate consecutively 
over fast degrees of freedom in the action (||) , which is de- 
fined with LOd being a high-energy cutoff. At each step of 
the RG procedure one has to eliminate fast modes W eie2 
in the N film, with either ft* > max(Z?<? 2 , e%— £2) > (for 
diffusons) or ft* > ma,x(Dq 2 , ei + e 2 ) > ^ (for Cooper- 
ons), and fast (with tt* > ui > fl) fluctuations of the 
order-parameter phase 9 on the SC island (where f2 is 
the running infrared RG cutoff). The above integration 
results in a correction to the action of slow variables 
proportional to the increment of the "logarithmic time" 
A£ = ln(Q* /Q). The structure of the boundary term in 
the action (^) is not reproduced under the RG Q, in- 
stead higher-order terms Tr(QsQ) n are generated, which 
are all relevant in the case of strong proximity effect. The 
full boundary action can be written in the form p2j 



-i7r 2 5 V 7 n(C)Tr(QsQ)". (11) 



Abound — ^ Sn 



At the energy scale md, the multicharge action (g_lj) re- 
duces to the last term in Eq. (^), that is, 71 (0) = Gt/^9 
and 7n>2(0) = 0. Integration over all fast modes down to 
the energy scale u> = u>de~^ defines the action (jll]) with 
the set of parameters 7 n (C)- 

In the model jl£| of large X n discussed above, separa- 
tion of scales was possible: at relatively short time scales 
the term S2 oc 72 Tr (QsQ) 2 was generated under the 
action of the RG [El, that lead to a constant value of 



Ga 1 (other 7„>3 were still small). At longer time 
scales fluctuations of 9 became important, being deter- 
mined by the action (^3j). In the full problem considered 
now, all parameters 7„ are important, and all types of 
fluctuations should be considered simultaneously. The 
corresponding RG equations for the "charges" 7 n (C) had 
been derived in Ref. Q for the case when the SC is- 
land is connected to an external circuit and its phase is 
fixed by an external bias V through the Josephson rela- 
tion d9/dt = 2eV. In the absence of an external contact, 
all phases do fluctuate and these fluctuations lead to an 
additional logarithmic renormalization of the action. 

To calculate the corresponding correction one has to 
average the action ([ll]) over high-frequency fluctuations 
of the phase 9 . Its propagator is determined by the Gaus- 
sian part of the same action and has the form: 



21K 



(12) 



where Ga(( u ) = 1 671-5 X^=i( -1 )" n2 72n(C0 is the scale- 
dependent Andreev conductance, £ w = hi(ujd/uj), and 



IE = 



2coth(w/2T) 1 



-1 







(13) 



Let us consider the n-th term of the action (|l l|) , S n — 
—i-K 2 g^ n r Tr{QsQ) n , and find the one-loop correction 
from it, A^7fc, to the parameters 7fc due to phase fluc- 
tuations with frequencies ui in the range fl < \lu\ < £1*. 
This correction will be proportional to AC = In(f2*/f2). 

To find it, one has to separate 9 into the sum of the slow 

part 9 (with frequencies < ft), and the fast part 9' 
(with < |w| < fi*), expand S n to the second order 

in 9' and average over it with the correlator ([l2|), taking 
into account all possible pairings. There are two types of 
pairing: two fast variables 9' can belong to the same ma- 
trix Qs [cf. Eq. (10)], or to two different Qs- The first 



case is simple: on averaging over 9' each of the terms 
S n reproduces its structure, and the whole effect of fluc- 
tuations is to modify 7n : A^7 n = — 2n7 n A£/7rG J 4(C)- 
The pairing of the second type involves averaging over 
the fast variables 9'(ti) and #'(£2) which belong to dif- 
ferent Qs- Therefore these two Qs subdivide the whole 
cyclic product (QsQ) n under the trace into two parts: 
Tr(QsQ)" = ■&Qs(ti)A(ti s t 2 )Qs(<a)B(*a,*i), where 
A = Q(QsQ) p , B = Q{Q s Q) n ~ p ~ 2 , and p is an inte- 
ger number which distinguishes different pairings. An 
example with n = 8 and p — 5 is shown in Fig. ^. For 
certainty we will assume that the fast energy lo runs over 
the matrix A (the opposite situation can be treated anal- 
ogously by replacing p «-> n — p — 2). Within the RG 
precision, fast bulk fields Q entering the product A can 
be substituted by A(w) = Kq{lj)t z (we will explain this 
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below) while slow matrices Qs entering A can be taken 
at coinciding times t m t\ m t<i. 




FIG. 1. An example of averaging (S n ) over phase fluctua- 
tions: n — 8, p — 5 (see text for details). Dots denote Qs(t), 
straight lines denote Q(t, t'), and the zigzag line stands for the 
correlator (^) of fast phases taken from Qs(ii) and Qs{t2)- 
The fast energy uj runs over the zigzag line and the part A 
of the diagram, the corresponding elements being shown by 
thick lines. 

Averaging is performed with the help of Eq. ([To]) and 
the following algebraic relations (which can be proved by 
direct calculation): 



6' ■ e** A (e* i8 A e ±i8 A )* ■& 
(e^Aoe^Ao)^^-^) 



0, 



(14) 
(15) 



According to these equations, only pairings with odd 
p = 2k — 1 produce a nonzero contribution. For them 
we have 

M k = (QsHQsAf^Qs) = (-l^^Qs, (16) 

where fast 0' are taken from the first and the last Qs 
under the correlator, and the sign (— l) k results from 
commutation of Pauli matrices in the Nambu space. 
Then, collecting all combinatorial factors, one gets AS* = 
(~l) k in~/ n (4irgA(/G A )TT(QsQ) n - 2k ■ We see that the 
paring with p — 2k — 1 modifies S n -2k with A^j n ^2k = 
(-l) fe+1 4n7„AC/7rGA(C)- Thus, under the action of the 
RG, the term S n generates the descending series of terms: 
S n , S n _4, . . . The RG equation for the coefficient 

7„ reads 



A 7n 



2AC 
'^GaJO 



2]T(-l) fe (n + 2fc) 7n _ 



2k 



k=l 



(17) 



Finally, we have to explain the validity of the substi- 
tution Q ^ A in the "fast" product A. Such reduction 
implies that no slow diffusons or Cooperons could be con- 
nected to the part A of the diagram on further steps of 



the RG. This statement is evident for Cooperons, but not 
for diffusons. Indeed, Cooperons depend on the fast fre- 
quency ei + £2 = 2w, whereas diffusons depend on e\ — £2 
which can be slow, and thus might contribute to the RG 
equation on further steps. Actually, however, they do 
not. The situation here is very similar to the one de- 
scribed during the derivation of the A term in the RG, 
cf. Sec. V of Ref. [^2): each diagram containing a dif- 
fuson in the "fast" product A produces several nonzero 
contributions, but their sum vanishes. 




FIG. 2. An example of diffuson pairing of a matrix Q from 
the fast part A of the diagram shown in Fig. |l] to the object 
Tr(QsA) 4 W . The double line denotes the diffuson propaga- 
tor. 

To prove this cancellation property, it is sufficient to 
consider the case when one Q from the part A is con- 
nected to some other object O by a diffuson, cf. Fig. ||. 
Since the pairing between S n and the object O is of 
diffusonic character, the latter ought to have the form 
O oc Tr (Qs-A) 2m W with even number of Qs- Pairing 
is performed with the help of the contraction rules from 
Ref. [^2| and yields two terms: 

Tr Qs(h)BQ s (t 2 ) [(A<5 s ) 2fc+2 " l_1 - (AQs) 2 ^ 2 ™- 1 ] A. 

(18) 



To obtain the final result, the expression (|18j) should be 
averaged over phase fluctuations. Note that Eq. ([l6]) is 
valid for any integer k. Therefore, averaging ( p^| ) with 
the help of Eq. (|16| ) yields A4+m — A4-m = 0, q. e. d. 

To simplify the RG equation, we introduce, following 
Ref. ]22] | , a function of an auxiliary continuous variable x 
according to the definition u(x, £) = J2n°=i n ln(() sinnx. 
Then the full RG equation for the function u(x, £) reads: 



+ uu x + A(C) w(§, C) siux 



TrG A {() 



(19) 



where !F[u(x 1 £)] = [w(x, £) tanx — u(§, C) secxj^, and the 
initial condition is it(x, 0) = (Gx/^g) sinx. The func- 
tional RG equation with the L.H.S. of Eq.(|i"9|) and zero 
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in the R.H.S. was derived in p2\ . In the present problem 
the R.H.S. describes the effect of phase fluctuations; it 
is obtained as the Fourier-transform of ([l?]). The scale- 
dependent subgap conductance Ga(C) is determined by 
the solution of Eq. @ as G A (() = 4tt#u x (§, ()• 

To find the parameter C(0) with exponential accu- 
racy, we integrate Eq. (|l^) together with Eq. (fit) for 
A(£). Written in the rescaled variables s = C,/2ir^fg : 
w(x,s) = 2n^/gu(x, £), and A = A/A s , Eq. ( |l9| ) acquires 
the form 



ww x 



^ Wx(f,s) 



A(s)w(f,s) since, (20) 



with the initial condition w(x,0) - 
G T /2^/g > 1. The solution of Eq. 



Asmx, where A = 
pp|) weakly depends 



on the ratio A(0) = X n /X g which is assumed to be not 
very large. At s <C 1, the function w(x, s) is close to the 
solution of Eq. ( ppf ) with zero R.H.S., which, at s > A^ 1 , 
is given by w(x, s) » x/s for x <E (0, n). As s grows, the 
R.H.S. terms become increasingly important, and even- 
tually reduce Ga(() — 2y/9 w x{^,s) down to the value of 
the order of 1 at the critical value of £* = 2-K^fgs c . The 

value of s c ~ 1 was determined, for several values of A(0), 
via numerical solution of Eq. (120) in the limit A — > oo: 



1.25, for A(0) = 0; 
1.17, for A(0) = 1; 
1.13, for A(0) = 2. 

As it was explained above, the value of (* 
determines C(0) with exponential accuracy: 

C(0) = 



(21) 
2ir^/gs c 

(22) 



where pre-exponential factor B is left undetermined (to 
find it, the two-loop RG calculation would be neces- 
sary). Combining Eqs. (g), @ and (p2|), we find that 
macroscopic superconductivity in the array is destroyed 
when the distance b between neighbouring islands be- 
comes larger than the critical value b c (g) defined as 



b c (g) 



2 V 



1/2 



(23) 



Finally, inverting relation (||^) one finds the result (0) for 
the critical film conductance g c (b), with 



0.64, for A(0) = 0; 
Gc = 4 = { 0.73, for A(0) = 1; 

0.79, for A(0) = 2. 



(24) 



&ndd = d-ir 3 / 2 13 1 / 2 g 1 / i /2s 1 c /2 . 



At <? > g c (6) macroscopic superconductive transition 
occurs at T > 0. Close to the critical point (jl]), at 
b < b c (g), the transition temperature is primarily de- 
termined by the temperature dependence of J(T) which 



is reduced compared to ,7(0) since couplings J(r) are 
strongly suppressed at r > Lt- At low temperatures 



J(T) = J(0) 



]n(L T /b) 
bx(L T /d) ' 



(25) 



The same expression, with Lt replaced by Ijj = 
■\/Trhc/eH, determines ,7(0, H) in the presence of a weak 
transverse magnetic field H. To find the critical temper- 
ature T c (b) and the T = value of the upper critical field 
H c2 , one uses Eq. (§) together with Eqs. (§|), ©■ The 
result is that both T c (g) and H C 2(g) scale in the same 
way, and drop fast at b — > b c (g): 



In- 



T* 



In 



$ 



21n(6/d) 



ff c2 6 2 bl{g)/b*-r 



(26) 



hD/b 2 , and $o = hc/2e is the flux quan- 



where T* 

turn. The form of Eqs. (26) suggests that the behavior 



of the array near the quantum critical point b = b c (g) 
can be understood in terms of the standard BCS-type 
theory with a macroscopic effective Cooper attraction 
A c ff, which vanishes at b — > b c (g). Eq. ( p6|) is valid for 

I/O 

b/b c (g) > (21n(6/<2)) / ~. This inequality ensures that 
T c is small compared both to T* (under this condition 
the proximity coupling is long-range) and to h/C(0). The 
latter condition allows one to approximate C(T) by C(0) 
while deriving Eq. (|2^). At shorter b -C b c (g)y/X^, the 
transition occurs at T c ~ T*. Here Lt c ~ the MFA is 
not applicable and the transition is governed by thermal 
fluctuations. Similarly, the characteristic scale of mag- 
netic field which affects superconductive state in such an 
array is just &o/b 2 : at higher fields formation of a super- 
conductive glass state is expected. Additional limitation 
for our results from the higher-T side is due to the ne- 
glect of quasiparticles inside SC islands, the correspond- 
ing temperature scale being T par i ty = Age / InivViAgc)- 
On the other hand, the conclusion about metallic na- 
ture of the system's state at b > b c (g) is limited to the 

2 

temperature scale T > T\ oc ~ ?kj t ;e _7r 9 , since we ne- 
glected weak localization effects (the value of g entering 
this estimate is defined at the length scale of the order 
of d) . Determination of weak- localization corrections to 
the conductance of an array with b > b c (g) would need 
the account of dephasing effects due to fluctuations of 
island's phases. This is an interesting problem which we 
left for the future studies. 

Finally, we comment briefly on the similar problem of 
small superconductive grains of radius d immersed into 
a 3D metal with bulk resistivity p. In this case the sim- 
ple method of Ref. |lq| can be used. The characteris- 
tic Coulomb energy Ti/C(0) « fru>d exp(— ttTl/ Ac 2 Ra) in 
the tunneling limit Rt 3> Rn', here Ra — R t /Rn and 
Jijv = p/fatd. Proximity coupling energy is estimated as 
J(r) ~ {h/e 2 R T ) 2 /r 3 is. The point of the SC-M transi- 
tion is then given by Tip/l6e 2 dR T = 31n6/d. 
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In conclusion, we developed a theory of quantum 
superconductive-metal transition in a 2D proximity- 
coupled array. The critical resistance Rn c is non- 
universal and small compared to the quantum resistance. 
Near the quantum critical point the system behaves as 
a BCS-like superconductor with the effective Cooper at- 
traction constant vanishing at R, n — > Rn c - 
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